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Abstract. 

We briefly examine the properties of dense plasmas characteristic of the 
atmospheres of neutron stars and of the interior of massive white dwarfs. These 
astrophysical bodies are natural laboratories to study respectively the problem of 
pressure ionization of hydrogen in a strong magnetic field and the crystallization of 
the quantum one-component-plasma at finite temperature. 
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1. Introduction 

The quest for an accurate description of the thermodynamics properties of dense plasmas 
has represented a thriving domain of research since the seminal review by Baus and 
Hansen (1980). The interiors of many astrophysical bodies are characterized by strongly 
correlated ionic and electronic plasmas, with respective classical and quantum coupling 
parameters Tj = fi(Z\e) 2 j a\ and r s = a- x j \a§Z\^) varying over several decades. Here, 
«o = h 2 /(m e e 2 ) denotes the electronic Bohr radius, a; = AtcV) 1 ^ is the mean 

inter-ionic distance, and (3 = 1/(/cb7 1 ). The correct description of the thermodynamic 
properties of these astrophysical bodies, which determine their mechanical and thermal 
properties, thus requires the knowledge of the equation of state (EOS) of such plasmas. 
In this short review, we focus on the case of neutron stars (NS) and massive white dwarfs 
(WD), which exhibit two particularly interesting problems in the statistical physics of 
dense matter. 

2. Ionization equilibrium of a hydrogen plasma in strong magnetic fields 

Most of neutron stars are characterized by magnetic fields B ~ 10 n -10 13 G, whereas 
some of them (so-called magnetars) are thought to have B ~ 10 14 -10 15 G. Although 
huge by the terrestrial standards, the magnetic energy of a NS, ~ R 3 B 2 /6, represents 
only a tiny fraction of its gravitational binding energy Eq ~ GM 2 /R 2 . 

The photospheric properties of a NS are characterized by temperatures T ~ 10 5 - 
10 7 K (depending on the age and mass of the star) and densities p ~ 10~ 2 -10 4 gcm~ 3 
(depending on T and B). The cooling rates of these stars are entirely determined by 
the relationship between the photospheric and the interior temperature profiles. The 
emitted spectra of these stars can be strongly affected by the presence of bound species in 
the photosphere. Since the atmosphere can be composed essentially of hydrogen accreted 
from either the interstellar medium, the supernova remnant or a close companion, 
the determination of the temperature profiles and spectra thus requires an accurate 
description of hydrogen ionization in a strong magnetic field. 

The quantum-mechanical properties of protons, free electrons and bound species 
(hydrogen atoms and molecules) are strongly modified by the field, which thereby affects 
the thermodynamic properties of the plasma. The properties of matter in a magnetic 
field under NS conditions have been reviewed recently by Ventura and Potekhin (2001) 
and Lai (2001). We refer the reader to these reviews for the detailed descriptions of 
these properties, and only the most recent results will be outlined in the present paper. 

The transverse motion of electrons in a magnetic field is quantized into Landau 
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levels. The energy of the nth Landau level of the electron (without the rest energy) is 
m e c 2 (v / l + 26n — 1), which becomes huj c n in the non-relativistic limit, where 

eB 

Hu c = h = 11.577512 keV, (1) 

m e c 

is the electron cyclotron energy, 

b = huj c /m e c 2 = B 12 /U.U (2) 

is the field strength in the relativistic units, and B 12 — Bj (10 12 G) is a typical magnetic- 
field scale for NS conditions. 

The atomic unit for the magnetic-field strength is set by Tiuu c = e 2 /ao, i.e. 
Bq = (m e c/he) x (e 2 /a ) = 2.35 x 10 9 G. It is convenient to define a dimensionless 
magnetic-field strength 

7 = B/B = b/a) , (3) 

where af is the fine structure constant. 

For 7 ^> 1, as encountered in NS conditions, the electron cyclotron energy is much 
larger than the typical Coulomb energy, so that the properties of interacting particles, 
protons, atoms, molecules, are strongly affected by the field. The ground-state atomic 
and molecular binding energies increase with B as ~ In 2 7. The atom in a strong 
magnetic field is compressed in the transverse directions to the radius ~ a m , where 

a m = (kc/eB) 1/2 = 7^ 1/2 a = 2.56 x KT 10 ^ 172 cm (4) 

is the quantum magnetic length, which becomes the natural length unit in the plasma 
instead of ao- 

The thermal motion of atoms causes the Stark effect due to the electric field 
induced in the comoving frame of reference. At 7 ^> 1, this effect strongly modifies the 
atomic properties: the atom becomes asymmetric, and its binding energy and oscillator 
strengths depend on the velocity (Potekhin 1994). Two classes of the atomic states 
arise: so-called centered and decentered states; for the latter ones the electron-proton 
separation is large and the binding energy relatively small (Vincke et al 1992; Potekhin 
1994). 

The formation of molecules is also strongly modified in a strong magnetic field. 
Because of the alignment of the electron spins antiparallel to the field, two atoms in their 
ground state (m = 0) do not bind together, because of the Pauli exclusion principle. 
One of the two H atoms has to be excited in the m — — 1 state to form the ground state 
of the H 2 molecule, which then forms by covalent bonding (e.g., Lai, 2001). 

As long as T < hu c /k B = 1.343 x 10 8 B 12 K and p < p B « 7.1 x 10 3 B% 2 gem" 3 , 
the electron cyclotron energy hw c exceeds both the thermal energy k B T and the electron 
Fermi energy A; B T F , so that the electrons are mostly in the Landau ground state, - i.e., 
the field is strongly quantizing. In this case, typical for the NS photospheres, the electron 
spins are aligned antiparallel to the field. 

Proton motion is also quantized by the magnetic field, but the corresponding 
cyclotron energy is smaller, hu cp = hu c m e /m p . 




Figure 1. Ionization isotherms at B = 10 G and two values of T. Solid lines: total 
fraction of atoms /h = nn/n ; dash-lines: fractions of ground-state atoms; dot-dashed 
lines: molecular fraction /h 2 = 2nn 2 /no, where no is the total number density of 
protons (free and bound). Dotted lines: zero-field case 



Quantum-mechanical calculations of the binding energies and wave functions of 
hydrogen atoms in any states of motion in the strong magnetic fields have been carried 
out only recently (Vincke et al 1992; Potekhin, 1994). Using these results, Potekhin, 
Chabrier and Shibanov (1999) derived a model which describes the thermodynamics of 
an interacting (H 2 , H, H + , e~) plasma in a strong magnetic field. This model is based 
on the framework of the free energy model developed by Saumon and Chabrier (1991, 
1992) for pressure ionization of hydrogen at B — 0, but generalizes it to the strong- £? 
case, taking into account the quantum-mechanical effects caused by the thermal motion 
of atoms across the magnetic field. 

Potekhin et al (1999) calculated the ionization equilibrium and EOS at 7x 10 11 G < 
B < 3 x 10 13 G. In this paper, we extend these calculations up to B = 10 14 G, typical for 
the magnetars. We calculate the number densities of atoms (tih) and molecules (nn 2 ) 
from the equations 

A p A e (2vra n 2 V2 



fin = n e -z — -3- 
A 3 



1-e 



-/3hoj c 



Z w e A , nn 2 =nl{\uV2fZ w2 /Zl, (5) 



where n e is the electron number density, \j = (27rf3h 2 /rrij) 1 / 2 is the thermal wavelength 
of the particle j (j = e, p, H), Z w and Z w2 are the internal partition functions for H 
and H 2 , respectively, and A is a correction factor, which takes into account electron 
degeneracy and filling of the excited Landau levels. The formulae for Z w , Z w2 and A 
are given in Potekhin et al (1999). 



Dense astrophysical plasmas 5 




log p (g cm" 3 ) log p (g cm" 3 ) 



Figure 2. EOS of partially ionized atomic hydrogen in the strong magnetic field (solid 
lines) compared with the EOS of a fully ionized ideal electron-proton plasma (dotted 
lines) and the EOS of a partially ionized hydrogen at B = (dashed lines). Left panel: 
B = 10 12 G, logT[K] = 6.5,6.0, and 5.5. Right panel: logT[K] = 6.5, B = 10 12 G 
and 10 14 G. The vertical lines correspond to the density above which excited Landau 
levels become populated 

Figure 1 displays the logarithm of the fraction of H and H 2 for a magnetic field 
B = 10 14 G, for two isotherms. The solid line represents the total fraction of atoms 
/h — nii/no (no = % + «)) + 2nn 2 ) in all quantum states, whereas the dashed line 
shows the fraction of atoms in the ground state. The dotted line displays the zero-field 
result. As seen in this figure, the strong magnetic field favours atomic and molecular 
recombination. Since the binding energies of atoms and molecules increase and Tp 
decreases in the strongly quantizing magnetic field, pressure ionization occurs at much 
larger densities than for the field-free case. 

Figure 2 displays the resulting EOS. The left panel displays three isotherms at 
B = 10 12 G, typical of ordinary NS. On the right panel, we compare the isotherm 
T = lO 6 5 K for B = 10 12 G and for the superstrong field B = 10 14 G. As expected from 
the previous discussion, the EOS in a strong magnetic field differs significantly from the 
field-free case in the region of partial ionization. At very high-density, excited Landau 
levels become populated, due to the increase of the Fermi energy, and the zero-field case 
is recovered. 
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3. Crystallization of white dwarf cores 

Massive white dwarfs (1.2 M < M < 1.4 M = M Ch , where M Ch denotes the 
Chandraskhar mass) are C 6+ /0 8+ plasmas with central density and temperatures 
Or < ^~' 10 8 -10 9 10 6 K. Under these conditions, the ion zero-point energy E oc 

Mlp, where VL P = (3Z?e 2 /M i af) 1/ ' 2 is the ion plasma frequency, exceeds the classical 
thermal energy k B T (Chabrier et al 1992), i.e. rj = %VL P /k B T 3> 1. Collective diffraction 
effects thus modify the classical Coulomb interaction. The melting values of the coupling 
parameters T m for the classical OCP and Rs, m , where R$ = a-J (h 2 / MiZfe 2 ) is the ion 
quantum plasma parameter, for the quantum jellium model at zero-temperature, have 
been firmly established: T m = 175 (Potekhin and Chabrier 2000), and Rs,m — 160 for 
bosons, Rs,m — 100 for fermions (Ceperley and Alder 1980). The melting curve at finite 
temperature, however, i.e., at n ^ and T^0, still remains poorly determined. 

A simplified determination of this curve is based on the Lindeman critical 
parameter interpolation between the zero-temperature and the classical melting values 
(Mochkovitch and Hansen 1979; Chabrier 1993). These calculations, however, are based 
on a harmonic description of the phonon mode spectrum and thus do not include non- 
harmonic effects. More recently, Jones and Ceperley (1996) performed Path Integral 
Monte Carlo (PIMC) simulations to try to determine the phase diagram more correctly. 
They found a maximum melting temperature almost a factor of 2 larger than the one 
determined by Chabrier (1993). 

In this paper, we present preliminary results based on similar PIMC calculations 
at finite temperature. The system under consideration consists of N identical, but 
distinguishable particles (Boltzmannions) with a mass M and charge Ze in a volume V 
at temperature T. 

The partition function of the quantum system is given by the trace of the N- 
body density matrix pn(R, R'', P) = (R\e~ l3HN \R!) . Here R denotes the 3iV coordinates 
{fi}i=i t ... t N- Using the Trotter formula, this partition function can be rewritten exactly 
in terms of the density matrix pn(R, R'] t) with r = f3/P (Ceperley 1995) as 

Qn,v,t(P) = J U.^ =1 dR a p N (R a , R a+1 ; r) (6) 
In the P — > oo limit, the exact density matrix 

e"^ = limp^e-^" x e~^] p (7) 
is recovered by using the approximate expression 

p N (R a , R p - r) = p° N (R a , Rp- t) exp [-^ [v c {R a ) + V C {R P )} 

where K N and Vjy denote the kinetic and potential parts of the Hamiltonian, 
respectively, p° N (R a , Rp;r) = (47rr)~ 3Ar//2 exp[— (R a — _R /3 ) 2 ]/(4r) is the density matrix 
of free particles, and V c is the classical potential energy of the system. 

In terms of the path integral formalism, the particle is defined by its trajectory 
in "imaginary time" Pt, through "polymers" composed of P monomers connected by 
"strings" of stiffness M-JTit. 



(8) 
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Figure 3. Internal energy of a quantum OCP at Rs — 200. Solid lines: PIMC 
calculations; left dash-line: Hansen- Vieillefosse (1975) fi, 4 -expansion; right dash line: 
analytical fit of the energy at large T (Douchin and Chabrier 2002). 



We have performed simulations along three isochores Rs = 1200, 350, 200 and 
five isotherms on each isochore around the estimated melting temperature. In order 
to estimate the finite-size effects, we used 16, 54, 128 and 250 particles. For each 
simulation we take rj/P =0.05-0.1. We then parametrized the fluid and solid internal 
energies, including the finite-size corrections. The melting temperature was estimated 
both from free energy comparison and from a dynamical criterion, namely the mean 
square displacement in unit of the nearest-neighbour distance ((Sr 2 ) /d 2 ) 1 ^ 2 . This value 
is finite in a solid and diverges in a fluid. 

Figure 3 displays the internal energy of the quantum fluid for Rs = 200. The 
solid line is the result of the PIMC simulation, with N = 54 particles. The dash line 
displays the Hansen- Vieillefosse (1975) fo 4 -expansion, which diverges when the thermal 
wavelength A* is of the order of the inter-ionic distance a;. At large T, i.e. large r] 
for a given Rs, the energy tends towards the zero-temperature limit (Chabrier, 1993). 
The dash-line displays an analytical fit of the PIMC energy at large F (Douchin and 
Chabrier, 2002). 

Our preliminary results yield melting temperatures lower than the ones obtained by 
Jones and Ceperley (1996), about a factor of 2 for the maximum melting temperature. 
This difference stems from a better correction of finite-size effects, which are known to 
stabilize the solid. First of all, we explored more in detail the N-dependence of our 
PIMC simulations by conducting calculations with a larger number of particles (up to 
N=256). Second of all, we found out that the dependence on F and Rs of the finite- 
size effects for the liquid phase is more complex than the simple one used by Jones 
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and Ceperley (1996). This is particularly important for the extrapolated values of the 
energy around the turning point. The obtained maximum melting temperature and 
corresponding density are (Douchin and Chabrier, 2002) 

T max = 8810 AZ 4 K; p m = 1280 A 4 Z 6 gem" 3 , i.e., R s « 235 . (9) 

These calculations show that, although anharmonic effects are non-negligible both 
in the solid and in the liquid phase, they are comparable in both phases and almost 
cancel out, so that the melting curve lies close to the one estimated from the harmonic 
spectrum. 

4. Conclusion 

In this short review, we have considered two different problems related to dense plasma 
physics as encountered under specific astrophysical conditions. Because of the presence 
of a strong magnetic field, the quantum internal levels of atoms and molecules are 
quantized in Landau orbits, and the field raises the binding energy of these species, 
favoring recombination over dissociation and ionization compared with the field-free 
case. This modifies the ionization equilibrium and EOS of the dense plasma. 

We also considered the crystallization of a quantum fluid of Boltzmannions at finite 
temperature with PIMC simulations, taking into account the correction due to finite- 
size effects. We found that the crystallization diagram lies close to the one based on an 
interpolation of the Lindeman criterium between the classical and the zero-temperature 
limits. These preliminary results need to be confirmed by further calculations. Such 
work is under progress. 

These two examples, and many other not mentioned in the present review, stress 
the need for detailed calculations of the properties of dense plasmas under extreme 
conditions for astrophysical applications. 
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